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Abstract 

This   report  examines   design   considerations  of  a  flow  system  for 
transient  electrical   laser  experiments.      After  a  brief  review  of  corona 
discharge   gas   flow  interaction  phenomena,    the   choice  of  the   size   and 
operating  conditions   of  the   discharge   cell   are   discussed.      Once   those 
parameters   are  established,    the   flow  system,   which   consists  essentially 
of  a  Ludwieg  tube,    is   designed.      Attention  is  paid  to   flow  starting 
times,    to   the  available   time  of  steady  state   flow,    and  to  the  effects 
of  the  boundary   layer  on   the   core   flow.      Suggested  experiments   and  in- 
strumentation  are  given,   as  well,    in  this   report. 


VII 


I.       INTRODUCTION 

The   effect  of  gasdynamic  means   for  stabilizing  corona  discharges 
in  atmospheric  air  has  been   shown      '    '        to  be   significantly  important. 
By  choosing   the   right  turbulence   spectrum  and  flow  speed,    the  electri- 
cal power  fed  into   a  diffuse   corona  discharge  has  been   improved  by  as 
much  as   500   times   above   the  maximum  no-flow  power.      The   influence  of 
turbulence,   which  is  more   substantial  at  the  higher  pressures    (-1   atm) , 
is   to  promote  better  mixing  and  homogeneity  of   the  gas.      Thus,   by  gas- 
dynamic  means,   molecular  gases  are   able   to   accept  higher  diffuse   dis- 
charge power,   making  such  discharge   attractive    for  electric  gas   laser 
operation  because  of   the  promise  of  attractive  output  power  capabili- 
ties. 

The   temporal  nonuniformities   of   the  gas   density  created  by  the 
turbulence  may  cause  phase   distortions   in  the    laser  output  beam,    and 
thus   deteriorate   the   laser  performance,    (due   to   the   effect  of  the  phase 
distortions  on   the   far  field  intensity  distribution) .      However,    turbu- 
lent  corona  discharge   studies   show  that   the   strongly  decaying   turbulence 

(2  ,4) 

mainly   affects    the   discharge   at  the   region   close   to  anode   tips 

where  the  electric   field  is   strong  and  the  plasma  is    formed;    this   is 
the   corona  region  of  length  L   ,    and  has   a  strong  nonuniform  electric 

cl 

field  E      (Fig.    1) .      The  external   zone   is   a  quasi-neutral   zone  which  has 
a  length  of  L,    L     >>  L    ,    and  an  electric  field  E     which  is   considera- 
bly lower  and  more   uniform  than  E    .      Thus,   by  exciting   the    turbulence 
at  the   anode   tips   region   and  letting  the   turbulence   decay  and  become 
very  weak  at  the  external  zone,   where   laser  action  occurs,    the   density 
nonuniformity  problem  can  be   reduced   to  minimum. 


The   separation  of   the   discharge   region   into   two  main   regions  with  space 
scales   L     and  L.     (where   L     <<   L.)    and  electric   fields  E&  and  Efa    (where 
E      >;>  E   )    appears  to  be   very  useful   for  efficient  laser  operation. 
Those   regions  may  be   controlled,   almost  independently,  by  changing 
electrode  geometry    (anode   tips   radii,   electrode   separation,   electrode 
material,   etc.)    and  flow  characteristics.      It  is  possible   to   choose   the 
average  electric   field  in   the  external   region   so   that  the  E/N  value 
will  be  in  such  a  range   that  most  of  the  electron     energy  will  go  to 
excite   the   upper  laser  level;    as   for  example,   in  aCO.   laser  containing 
a  mixture  of  He/N^/co^   in  a.  volumetric  ratio  of   3/2/1,  more  than  80%  of 
the  electron  energy  goes   into  the  co_  (001)    and  N    (v=l)    vibrational 
levels   for  the  value   of  E/N   -  2   x  10     °  vcm  '      .      In  this  respect, 

an  analogy  can  be  made  between   a  c.w.   electron-beam-sustained  electric 

discharge  and  a   c.w.    gasdynamic  stabilized  corona  discharge;   both  cases 

(7) 

operate  in  the  Ionizer/Sustainer  mode    ,  which  means  that  there  is  a 

separation  between  the  electric  field  that  produces  the  electrons  and 
the  sustained  discharge,  where  the  electrons  acquire  their  optimal 
temperature  Te  (Te  is  determined  by  E/N)  for  laser  action. 

In  this  report  we  propose  an  experimental  study  of  the  effects  of 
gasdynamic  means  (mainly  turbulence  and  flow  speed)  and  different  geo- 
metrical parameters  of  the  electrodes  on  corona  discharges  for  laser 
applications. 

The  experimental  system  proposed  is  a  versatile  inexpensive  set-up, 
in  which  the  following  can  be  studied. 

Fluid  dynamics  measurements : 

Turbulence  spectrum      Flow  rates 

Temperatures 
Velocity  profiles        Pressures 


Electrical  measurements: 

V-I  characteristics 

probes  for  E,  Ne,  and  T  measurements 
Optical  measurements : 

Photography 

Schlieren  photography 

Interferograms 

Light  emission  spectroscopy 
Laser  performance  measurements  : 

Small  signal  gain 

Output  power  measurements 
Parameters  to  be  studied: 

Turbulence  intensity  and  spectrum 

Flow  rates  (speed  and  pressure) 

Gas  mixtures 

Discharge  electrical  energy 

Distance  between  electrodes 

Electrodes  geometry 

Preionizing  devices,  like  u.v.  radiation 

To  generate  the  flow  through  the  corona  discharge,  the  application 

(8) 
of  the  Ludwieg  tube    principle  has  been  considered  to  be  an  attractive 

proposition  because  of  its  simplicity,  low  cost  and  also  high  flow 

quality.   The  useful  flow  time  of  this  device  depends  mainly  on  the 

charge  tube  length  and  on  the  speed  of  sound  of  the  gas  in  the  tube  (see 

,  .  .  r-  msec 

Section  III).      For  air  this   time   is  approximately  5  - — ; —  .      For  a 

rc  J        meter   length 

charge   tube  of  the  order  of  4m  in  length,    the  net  useful   flow   time  will 


be   15-20  msec,   which  corresponds  to   -30   cavity   flow  times    (for  5   cm 
cavity  length  and  gas   speed  of  100  m/sec)  .      This  procedure  effectively 
simulates   cw  operation  without  requiring  the  heat  sink  electrodes/ 
large   compressors,   heat  exchanger  etc.,    that  would  be  necessary  for 

true   c.w.    operation.      All  other  relevant  characteristic   times   are  much 

(2) 

shorter  than   the   flow  time 

Section  II   deals  with  the  basic  design   and  sizing  of  the   test  cell. 
Once  the   cell  dimensions,    the  pressure  of  the  gas   and  its   flow  speed  are 
determined,    the   Ludwieg  tube   flow  system  can  be   designed.      The   detailed 
considerations  of  the  design  are  given  in  Section  III.      The  overall  ex- 
perimental set-up  is   shown  schematically  in  Fig.    2. 


II.       TEST    CELL   DESIGN 

The   test  section  design  is  basically  the  same   as   that  of  the  cell 

(2    3) 
used  previously  in  the  gasdynamic-discharge   interaction   studies 

Some  changes  have   to  be  made   due  to   the  use  of  this   cell   for   laser  ex- 
periments and  as  being  a  part  of  a  different  flow  system. 

1.      Cell  Operating  Conditions 

Gas  mixtures     Different  gas  mixtures   are  going  to  be  used  in 

this  system:      N   ,   He,    CO    ,   CO    ...   etc.      As   far  as   a  flow  system 

design  is   concerned,    the  difference  among   those  mixtures   is   their 

specific  heat   ratio  y   stnd  molecular  weight, 

will   consider  two  values   for  y : 

Y  =   1.4    (diatomic  molecular  gas) 

Y  =   1.67    (monatomic  gas) 

Gas  pressure     The  cell  will  be  operated  at  pressures   around  500 
torr  which  is  the  pressure   at  which  maximum  electrical  power   can  be 
dumped  into   the  corona  discharge   in  no   flowing  air   wi-fch  4cm     elec- 
trode gap        .      For  different  gases,    different  electrode  gaps   and  for 
flowing  gases   the  optimum  operating  pressure      varies,    thus   the  sys- 
tem will  be  designed  to  operate  in  the  pressure   range  of  200   <  p   « 
800  Torr    (lower  and  higher  pressures   can  be  achieved  as  well) . 

Flow   speed     The   gasdynamics-discharge    interaction  becomes   effec- 
tive at  flow  speeds  of  ~  100  m/sec        .      The  faster  the  flow   is,    the  more 
discharge  performance  improves.      The  cell  will  be  designed  to  have   a 
flow  of   about  Mach  number  M  =  0.5. 


2.  Cell  dimensions 

Cross  section  dimensions: 
Cell  height:    d  =  2.22",  5.5  cm  . 
Cell  width:    b  =  6.5"  ,  16   cm  . 
Dimension  along  the  flow  direction: 
Cell  length:    c  =  4"    ,10   cm  , 

Dimensions  d  and  c  are  the  same  as  those  in  Ref.  3;  however, 
the  cell  width,  b,  is  bigger  since  along  this  direction  the  small 

signal  gain  (s.s.g.)  is  measured  and  laser  cavity  will  be  aligned. 

( 6  9 ) 
The  s.s.g.  reported  in  the  literature   '    for  electrical  lasers, 

range  from  1%  to  4%  cm   ,•  16  cm  cavity  length  should  be  adequate  for 

gain  and  power  measurements. 

3.  Nozzle  throat 

The  cross  section  area  of  the  discharge  cell  is  S  =  b  x  d  = 

2 

88  cm  .   In  order  to  obtain  M  =  0.5  in  the  test  section,  the  cross 

section  area  A*,  at  the  throat  for  the  different  values  of  y,    are 
given  in  Table  I: 


Y    =    1.4 

Y    =    1.67 

A/A* 

1.34 

1.31 

A* (cm  ) 

66 

67 

h*(cm) 

4.1 

4.2 

Table  I 


Table  I  give  also  the  values  of  throat  height  h*,  for  a  one 
dimensional  nozzle. 

Turbulence  generating  plates  have  to  be  inserted  ahead  of 
the  electrode   tips.   Since  the  effective  cross  section  of  those 


plates  is  smaller  than  that  of  the  cell,  their  effective  cross 
section  should  exceed  the  values  of  A*  given  in  Table  I;  other- 
wise, the  flow  will  choke  at  the  plate,  prior  to  reaching  the  test 
section. 

The  turbulent  plate  used  in  previous  experiments  had  an 
area  blockage  of  the  order  of  50%,  which  is  much  above  the  34% 
allowed  for  this  design.  A  new  design  and  location  of  those 
plates  have  to  be  considered.   One  way  of  improving  the  plate 
blockage  is  to  mount  the  turbulence  generators  closer  to  the 
electrode  tips.   This  will  enable  a  weaker  turbulence  generation 
(thus,  a  smaller  blockage)  since  the  decay  of  the  turbulence  at 
the  electrode  tips  will  be  insignificantly  small.   One  way  of 
doing  it  is  shown  in  Fig.  3. 
4.   Electrical  power  heat  addition 

The  effect  of  adding  heat  to  a  compressible  flow  is  to  send 
compression  waves  into  the  flow  which  change  flow  parameters ;  heat 

addition  to  subsonic  flow  will  accelerate  the  flow  and  for  a  cri- 

(10) 
tical  value  of  heat,  q   ,  the  flow  will  be  choked. 

^cr 

o   -   qcr   =     1       (1  +  ym2)2  _ 
cr   C  T  y  -  '   2 

p  stag   2(y  +  1)M   1  +  — - —  M 


Typical  mass  flow  rates  through  the  cell: 

P 


m  =   pAV  = 


(R/M)T 


Av 


i2.  4 

2      x  88  x  2    x   10  ,    9     . 

for  N_:  m  =  = =    1  Kg/sec 

(8-32g  10   )«300 


i26 

2      x  88  x  2   x   10  rt   ,  .   _     . 

for  H    :  m  =  = -   0.14  Kg/sec 

8                                ,8.3   x  10   t    „„_ 
( : )*  300 


Typical  critical  heat  addition   values  are  given   in  Table   II 


2Cr 

^■cr    kg/sec 

q   ,  kw.For  typical 
cr    '  .    , 

m  values 

for  N2 

0.45 

140 

140 

for  H 

e 

0.4' 

620 

87 

Table   II 


The  power  supply  available  has  ~ 1  kw  power  output  which  is 
very  low  power  compared  to  q 


III.   LUEWIEG  TUBE 
1.   Introduction 

The  Ludwieg  tube  flow  system  is  shown  schematically  in  Figure  2. 
The  Ludwieg  tube  consists  of  a  long  tube,   charge  tube,  where  the  gas 
is  initially  compressed  to  the  desired  charge  pressure,  P  .   A  tran- 
sition section  channels  the  flow  from  the  circular  charge  tube  to 
the  rectangular  test  section  which  is  a  part  of  a  converging-diverging 
supersonic  nozzle.   The  gas  expanding  through  the  nozzle  is  discharged 
through  a  short  diffuser  into  a  dump  tank.   A  diaphragm  separates  the 
driver  and  driven  sections.   The  diaphragm  can  be  located  either  down- 
stream or  upstream  of  the  nozzle. 

The  flow  is  initiated  by  the  rupture  of  the  diaphragm.   In  the 
case  of  an  upstream  diaphragm,  a  backwards  facing  centered  expansion 
fan  propagates  into  the  charge  tube.   This  fan  is  reflected  at  the 
closed  end  of  the  tube  and  returns  back  to  the  nozzle.   This  is  called 
the  first  flow  cycle;  a  simplified  wave  diagram  is  shown  in  Fig.  4(a) . 
Between  the  time  when  the  tail  of  the  expansion  wave  leaves  the  nozzle 
entrance  and  when  the  head  of  the  reflected  fan  reaches  the  nozzle 
again,  the  stagnation  conditions  ahead  of  the  nozzle  remain  constant  and 
a  clean  steady  flow  may  be  established  in  the  nozzle.   A  second  flow 
cycle  (and  others)  follows  the  first  one  but  the  stagnation  parameters 
are  different.   The  duration  of  the  steady  flow,  for  a  certain  gas,  is 
primarily  a  function  of  the  charge  tube  length  and  we  are  interested  to 

make  it  as  long  as  possible;  there  is  a  limit  on  the  length  of  the  tube 

(11  12  13) 
set  by  the  growth  of  the  boundary  layer  "     '    .   The  boundary- layer 

growth  sets  also  a  limit  on  the  diameter  of  the  tube  and  on  the  flow  Mach 


number,  M  ,  in  the  charge  tube.   Subsequent  to  the  diaphragm  rupture, 
a  shock  wave  and  a  contact  surface  proceed  downstream  through  the 
driven  section.   Additional  waves  follow  the  contact  surface  to  adjust 
the  nozzle  test  section  pressure  to  the  freestream  pressure  correspondr 

ing  to  the  designed  Mach  number  expected  at  the  nozzle  exit  for  a  given 

(14) 
nozzle    .   The  steady  supersonic  flow  in  the  nozzle  is  not  established 

until  all  those  waves  are  swept  out  the  nozzle.   The  time  required  to 
get  the  flow  started,  that  is,  the  time  required- to  achieve  a  steady 
flow  from  rest,  in  the  supersonic  nozzle,  is  a  very  important  performance 
criterion.   It  is  important  to  minimize  the  start  time  in  order  to  maxi- 
mize the  duration  of  the  steady  state  conditions  for  a  fixed  charge  tube 

..  (14,15) 
length 

For  downstream  diaphragm  location,  the  wave  diagram  is  different 
from  the  case  where  the  diaphragm  is  located  upstream,  a  simplified  x  -  t 
diagram  is  shown  in  Fig.  4(b) .   As  soon  as  the  local  sonic  velocity  is 
reached  at  the  nozzle  throat,  no  further  expansion  waves  can  pass  through 
into  the  tube  and  an  expansion  fan  of  finite  width  passes  down  the  tube. 

For  an  upstream  diaphragm  location,  the  starting  times  at  the 

nozzle  exit  are  reported  to  be  shorter  than  those  obtained  with  the 

(15) 
downstream  diaphragm  location    .   However,  the  disadvantage  of  the  up- 
stream location  diaphragm  is  that  when  breaking  the  diaphragm, particles 
might  tear  off  and  disturb  the  flow  in  the  test  section. 
Charge  tube  diameter  and  length 

Economic  considerations  dictate  that  a  minimum  charge  tube  diameter 
be  used  for  a  given  test  section  size.   This  results  in  a  relatively 
small  contraction  (small  ratio  between  nozzle  throat  area  to  charge  tube 


10 


area)  which  requires  high  Mach  number,  M  ,  at  the  charge  tube.   This 
will  cause  the  boundary  layer  to  have  a  dominant  effect  on  nozzle  flow. 
It  is  suggested  to  have  the  I.D.  of  the  charge  tube  to  be  the  same 
as  the  dimension  b  of  the  test  cell,  i.e.  16  cm.   In  this  case  A./A* 
and  M^  are  given  in  Table  III. 


Y    =    1.4 

Y    =    1.67 

A4/A* 

3.05 

3.00 

M3 

0.2 

0.19 

Table  III 


If  a  tube  with  a  larger  diameter  is  available  (like  20cm) ,  it  will  im- 
prove the  performance  of  the  system,  since  M_  will  be  smaller.   It  is 
not  recommended  to  use  a  plenum  chamber  upstream  of  the  nozzle  inlet 
(to  slow  and  mix  the  core  and  boundary  layer  flows  emerging  from  the 
pipe)  since  it  will  be  shown  that  for  M_  =  0.2  the  disturbances  intro- 
duced by  the  boundary  layer  growth  are  small  while  the  starting  time  is 
going  to  be  much  longer  by  introducing  the  plenum  chamber. 

6*/d  is  given  as  a  function  of  ins  tan  aneous  aspect  ratio  of  the 


flow  s/d,  P  d  and  M  ,  in  the  form  of  5*/d 

(16) 


(P4d) 


(s/d) 


4/5 


by  Russel  ,   based 


on  Becker's  model   for  5*/d  ,   where  5*  is   the  boundary   layer  displace- 

ment  thickness,    at   a  distance   s   from  the  expansion  wave   head,    and  d  is 
the   charge   tube   diameter. 


for  M3  „  0.2     z  7-774/5 


=  5  x   10"4    (atm  -  m) 1/5 


(s/d) 

At  the  end  of   the   test  time,    s    =  2L  and  the   thickness  of  the  boun- 
dary layer  is   6    =    106*.      For  P.   =   1   atm,    L  =   3  m  and  d  =  0.16  m, 
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we  get  6*/d  -  1.3%,  6  -  0.02  m  and  5/d  -13%. 

The  perturbations  in  the  core  flow  due  to  the  boundary  layer, 

when  the  test  period  is  over,  are  (for  the  same  parameters  used  for 

6* 

— —  calculations) : 

d 

_AP         (P   d)1/5  AP 

—    — - =   4.3  x   10   *    ^   -  — -=•  =   6.5    x  10 

p  P 

3       (L/d)4/5 

-v —475   =   3'2  x  10  v3 

3    (L/d)4/i 


AT  Y   "    1   AP 

— -  = -r-  =  -  r~  x  6.5  x  10    =  -1.7  x  10 

T3  y  P3  1-4 

which  are   relatively  very  small  disturbances. 

The  effects  of  the  above  described  core  perturbations  on   the  down- 
stream flow  were  discussed  by  Russell  .      He  evaluated  the  boundary 
layer  along  the  walls  of  a  set  of  specified  nozzles.      The   results   for 
different  values  of  6*/d,    M     and  M    (Mach  number  in  the  nozzle)    are 


presented  in  Fig.    7  of  his  paper;  he  presents  the  effect 


lve   area  - 


(A/A*)eff  -   A/A*  _  .  .    »,, 

'a/a* as   a   functlon  of    °*/d  for  different  values  of  M     and 

M  where    (A/A*)  is   the   test  section  -   to  -   throat  area  ratio  with 

•  eri 

boundary-layer  correction. 

A/A*  is   the  geometric  area  ratio,   with  no   corrections. 
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For  the  value  of  P  d  =  0o5  atm-m,  Russell  found  that  with 

-4   x 
10   <__  o*/d  <_   1% ,  the  effective  area  change  parameter  has  a  finite 

value  determined  by  the  normal  steady  boundary  layer  correction  of 

the  nozzle  itself  (in  our  case  P  d  -  2  x  0.16  z.   0.3.   There  is  a 

o 

weak  dependence  on  pressure   level    (P.d)  )  . 

To   summarize: 

for  M  ■  0.5  Mach  number  flow  in  the  test  section,    d=16cm,    L=3m 

and  m     =  0.2 

<S*/d  =  1.3%  t  5/d  =   13% 

AP 

3  '  "3 

(A/A*)  ~   A/A* 

eff  m  (contribution   from   charge  tube 


A/A*  boundary   layer) 

3.      Test   time 

The  test  time,    (the  time   taken   for   the  head  of  the  reflected 
expansion  wave   to  reach  back   the  throat)    for  small  values  of  M     is 
given  approximately  by  2L/a   ,   where   a.    is   the  sound  speed  of  the 
undisturbed  gas   in  the  charge   tube. 

The  time  available   falls  only  slightly  with   increasing  the  Mach 
number  M_ 

Y  +  l 


i^{l  +  .L^M}2(Y 
o        1  +  M_  2  3 


t/t_   =   ;     ,    „        (1  +  -L^-^  M,}' 


where  t  =  L/a..   For  M_  =  0.2  and  L  =  3m,  the  calculated  values  of 

O  4  J 

t     and  t  are  given  in  Table   IV. 
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Y    =    1.4 

Y    =    1.67 

a    (m/sec) 

350 

1000 

t     (msec) 
o 

8.55 

3 

**o 

1.88 

1.9 

t (msec) 

16 

5.7 

Table  IV 


Starting  times 

The  starting  time  of  supersonic  flow  at  a  given  location  in  the 
flow  system  is  defined  as  the  time  elapsed  between  the  arrival  of  the 
head  of  the  expansion  fan  (downstream  diaphragm) ,  or  shock  wave 
(upstream  diaphragm)  and  the  instance  where  the  pressure  becomes  time  in- 
dependent, i.e.,  when  steady  flow  is  established.   For  an  upstream 
diaphragm  location,  the  starting  times  at  the  nozzle  exit  are  reported 

to  be  shorter  than  those  obtained  with  the  downstream  diaphragm  loca- 

<■   (15) 

tion 

a.   Downstream  diaphragm  location : 

According  to  ref.  15,  (they  used  a  Ludwieg  Tube  with  a  charge  tube 
tube  of  670  cm  in  length  and  I.D.  of  13.4  cm)  the  starting  times  at  the 
nozzle  exit  are  independent  of  the  initial  downstream  pressure  p  , 
except  in  the  range  above  performance  limit.   Starting  times  increase 
with  diffuser  length.   The  starting  times  upstream  of  the  nozzle  throat 
do  not  depend  noticeably  on  diffuser  length,  and  are  of  the  order  of 
0.5  msec,  where  the  starting  times  in  the  diverging  part  of  the  nozzle 
are  of  the  order  of  6  msec.   The  experimental  results  are  presented  in 
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dimensionless  starting  time  t  /x  as  a  function  of  dimensionless  axial 


distance  x/d,  where  t  is  defined  as  &/a-v;  I'  -   /R*h*,R*  is  the  radius 
of  curvature  at  the  nozzle  throat,  h*  the  nozzle  height  and  a*  the 

speed  of  sound  at  the  throat.   At  the  subsonic  nozzle  section  t  /t-  =  3  anc^ 

c  s   f 

at  the  nozzle  exit  t  /t   ~  40.   In  our  case,  the  starting  time  in  the 
subsonic  section  will  be  t  "  1  msec  and  in  the  supersonic  section 
t   ~  8  msec,  for  N  .   For  other  gases,  the  starting  time  scales  in- 
versely to  a*, 
b.   Upstream  diaphragm  location 

The  starting  times,  at  the  exit  of  the  nozzle,  increase  with 
initial  downstream  pressure,  P  ,  and  increase  also  with  diffuser  length. 
The  starting  times  are  of  the  order  of  2  msec  at  nozzle  exit  and  ^4/^1 
should  be  at  least  of  the  order  of  10  in  order  to  start  the  flow 
Summary : 

It  is  recommended  to  use  the  downstream  diaphragm  location, 
mainly  because  the  starting  time  in  the  subsonic  portion  of  the  nozzle 
is  very  short  (  -  1  msec)  and  the  troubles  of  having  diaphragm  particles 
in  the  flow  are  eliminated.   For  supersonic  experiments,  it  might  be  an 
advantage  to  have  the  diaphragm  located  upstream.   The  diaphragm  holder 
unit  can  be  designed  so  that  there  will  be  a  possibility  to  mount  it 
either  downstream  or  upstream. 
Charge  tube  pressures  and  temperatures 

The  pressures  and  temperatures ,  P  and  T  ,  in  the  charge  tube ,  prior 
to  diaphragm  rupture,  as  a  function  of  pressure  P  and  temperature  T  , 
respectively,  where  the  iMach  number  of  the  flow  is  M  ,  is  given  by    : 
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Y  -  1   2 

T_    T  ^  T_      T   T    T      1  + M 

2     st  2   _  _3   st.   2   _  f £ 

T „  ~  T.   T  A    =  T.  T_   T  .    M    Y  -  1..  .2 
4    4    st     4   3     st    (1  + 


M     .  +  ^— ■; M„ 

2    3         2     2 


(1) 


Y 


-a-  i-V"1 

P4    T4 


(2) 


For  Y  ■  1.4,  M  ■  0.2  ■*  —  »  0.93 


i  +  X^A^ 


For  Y  =  1-67,  M  =  0.2  *  -2-  =  0.95  X 

The  corresponding  values  of  p4/p2  can  be  calculated  from  Eq.  (2)  . 

P.  as  a  function  of  P_  for  different  values  of  the  M_  is  shown 
4  2  2 

in  Fig.  5  and  Fig.  6  for  y  =  1*4  and  y  -   1.67  respectively.   The 
values  of  T  ,  for  y  =  1.4  and  y  =  1.67  and  different  values  of  M  are 
given  in  Table  V  assuming  T  =  300  k: 


M2 

0 

0.5 

1 

1.5 

2 

Y  =  1.4 

280 

266 

233 

193 

155 

Y  =  1.67 

285 

263 

213 

162 

122 

Table  V 


The  pressure,  P.,  in  the  part  of  the  system  downstream  the  diaphragm 
should  be  of  the  order  of  P  /100  (to  get  best  start  times  without  a 
diffuser) ;  i.e.  P   =  1  torr  (see  Sec.  5) . 
5,  Diffuaar,  dump  tank  and  valve 

The  starting  times,  as  discussed  before,  become  longer  with  longer 
diffusers,  both  for  downstream  and  upstream  diaphragm  locations. 
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For  a  fixed  charge   tube  pressure,   and  back  pressure  P..    being  in- 
creased gradually,    the   flow  at   the  diffuser  end  ceases  to  start  for  a 
certain  value  of  P..    while   further  increasing  P..  ,    the   supersonic   flow  at 
the  nozzle  exit  barely  starts.      Beyond  these   critical   conditions,    the 
flow  at  the  nozzle  exit  becomes   unstable.      The   limiting  experimental 
downstream  pressure   radios,   P-j/Pw   were  measured  by  Merkli  and  Abuaf   '       , 
both   for  upstream  and  downstream  diaphragm  locations,    as  a  function  of 
diffuser  length  L/D. 

It  was   found  that    (P. /P.)         for  diffuser  starting  is   -0.11   inde- 

1  4  cr  ^ 

pendent  on  L/D. 

(P.. /P.)    for  nozzle  supersonic  flow  increases  with  L/D  and  reaches 
1  4  cr  r 

a  maximum  at  about  L/D  =  6  where  (P, /P.)    =  0.22. 

1     4   cr 

That  means   that  within   the  available    time,    the  pressure   in   the  dump 
tank  should  not  increase   above   0.1  x  P.;    this   limiting  pressure   defines 
the  volume  of  the   dump  tank. 

The   rate  of  pressure   rise   inside   the   tank:      The   continuity  dictates 

apd 
that  in  the  dump   tank  V.  r——  +  m,   =0  where  V,  is   the  tank's  volume,    p. 


is   the  density  and  m,  is   the  total  mass   flow  into   the   tank.      By  differen- 
tiating the  equation  of  state,    P     =  p,  R  T,,   with  respect  to   time   and 
using   the   continuity  equation,   we  get 

3Pd         „     3pd       RTd   . 

—  '  RTd  at"  =  vT  md 

a 

P3   A3M3a3 
md  =  m3.  =  P3  A3   V3   =  R?J~ 


V3A3M3a3 

^d" T^ T 

3d 
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Where  T  is  the  available  testing  time.   Since  initially  P  =  P.  << 

0.1P,  we  may  write  AP,  =  P.;  from  the  last  equation: 
3  da 

n  7  2  —3 

V-     =    ( — )  (— )A„M„a^    <    2    x   10  x    nx8  xO.2    x   33000   x  20    x   10 
d  T3      Pd      3    3°3     - 

VJ   -   5.3   xlO      cm 
d 

It  is  suggested  to  install  a  valve  \zpstream  the  diapragm,  as 
shown  in  Fig.  3.   After  each  run,  before  the  ruptured  diaphragm  is 
taken  off,  this  valve  will  be  closed  and  the  gas  in  the  system  will 
be  saved. 

A  more  efficient  way  of  saving  gas  will  be  to  use  a  compressor 
to  compress  the  gas  in  the  dump  tank  and  recirculate  it  back  to  the 
charge  tube.   The  compressor  should  be  vacuum  tight.   A  simple  and 
inexpensive  compressor  is  suggested  in  Fig.  7. 
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IV.       INSTRUMENTATION 
Pressure   gauges 

0-5   atm.    range   for  measuring  P       . 
0-100  torr 

0-760   torr    y     Pressure  gauges   for  gas  mixtures  preparation 
0-2   atm 
0-100  torr   for  measuring  P. 

For  measuring  P_    and  P     during  the  experiment   itself,   after 
diaphragm  is  broken,   we  need       piezoelectric  pressure   transducers  with 
time   response  of  the  order  of  tens  microseconds,    in   the  range  of 
0-2   atm. 
Velocity  profiles^  and  turbulence   spectrum 

For  transient  flow   speed  and  turbulence  measurements,    a   laser- 
doppler  anemometer  or  a  hot-wire  anemometer  can  be   used.     The  avail- 
able  steady   flow   time,    10-20  msec,   might  be   too   short   for  appropriate 
averaging  the   turbulence   spectrum,   mainly  in   the   low   frequency   regime. 
Therefore   it  will  be  necessary   to  run  the   flow  in  a  continuous  mode    (with 
compressed  air)    for  turbulence   spectrum  measurements.      The   transient   flow 
fluctuations  may  not  be  averageable  but  their  ability  to   stabilize   the 
discharge  may  surpass   that  of  steady   flow  in  which   case  other  measurement 
forms   should  be   sought. 
Voltage    current  measurements;    probes 

The  voltage  and  current  have   to  be  recorded  during  the   experiment. 
The  simplest  way  of  doing   it  would  be  by  using  an  oscilloscope. 

The  probes   used  for  measuring  plasma  parameters  must  be  provided 
with  automatic  bias   sweep,    so   that  the   characteristics   can  be 
accomplished  within  the   flow  time. 

19 


Lasing  diagnostics 

A  probe  laser  for  s.s.g.  measurements. 

A  fast  response  detector  (microseconds  range) 

Power  meter. 
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DIAPHRAGM 


(3)    FIRST  FLOW  CYCLE 


EXPANSION  FAN 


X 


SHOCK   WAVE 
L  CONTRACT  SURFACE 

UPSTREAM  FACING   WAVE 


Fig.  4(a)     APPROXIMATE   WAVE  DIAGRAM  FOR   THE  FLOW  IN   A 
LUDWIEG  TUBE.    —   UPSTREAM  DIAPHRAGM  LOCA- 
TION. 
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DIAPHRAGM 


EXPANSION  FAN 


SHOCK   WAVE 
CONTACT   SURFACE 


Fig.    4(b)    APPROXIMATE   WAVE   DIAGRAM  FOR  THE   FLOW   IN   A 

LUDWIEG   TUBE.- DOWNSTREAM  DIAPHRAGM  LOCATION 
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Fig.    5    CHARGE  TUBE    PRESSURE,  P       AS   A  FUNCTION  OF  CELL 
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PRESSURE,  P       FOR  DIFFERENT   VALUES  OF   M         y  =1.4 
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P    (torr) 
4 


1800  r- 


1600  - 


1400  - 


1200  - 


1000  - 


800  - 


600  - 


400  - 


200  - 


X=1.67 


100  200  300  400  500  600 


700 


P2  (torr) 
Fig.   6    CHARGE   TUBE   PRESSURE,  P        AS   A  FUNCTION  OF  CELL  PRESSURE, 
FOR  DIFFERENT   VALUES  OF  M       y  =1.67 
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